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Abstract 



We report the first measurement of a structure dependent component in 
the decay K~^ ^^'^v^^y. Using the kinematic region where the muon kinetic 
energy is greater than 137 MeV and the photon energy is greater than 90 
MeV, we find that the absolute value of the sum of the vector and axial- 
vector form factors is \Fy + Fa\ = 0.165 it 0.007 it 0.011. This corresponds to 
a branching ratio of Si?(SD+) = (1.33 ± 0.12 ± 0.18) x 10"^ We also set the 
limit -0.04 < Fy - Fa < 0.24 at 90% c.l. 
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Typeset using REVT^ 



The decay K^^^'^v^'-y [K^^^) can proceed via two distinct mechanisms. The first, 
internal bremsstrahlung (IB), is a radiative version of the familiar ^ ^~^i>^ {^^2) decay: 
its Feynman diagram has a photon emitted from the external kaon or muon line. The 
second, structure dependent radiative decay (SD), involves the emission of a photon from 
intermediate states. SD is sensitive to the electroweak structure of the kaon and has been 
the subject of an extensive theoretical literature In recent years most of this has been 
in the framework of Chiral Perturbation Theory (ChPT) ||^. The differential rate in the 
rest frame can be written in terms of x = |^ and y = ^^^jj^^, where is the photon 
energy, i?^ is the muon kinetic energy, is the fi~^ mass, and is the mass: 
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In these formulas, Fy is the vector form factor. Fa is the axial form factor 0], a is the fine 
structure constant (1/137.036), Fk is the decay constant (159.8 ± 1.4 ± 0.4 MeV), and 
Tk^2 is the width of the K^2 decay. 

SD"*" and SD~ refer to different photon polarizations and these components do not mu- 
tually interfere. Both SD^ and SD^ can interfere with IB, however, resulting in the terms 
labeled INT"*" and INT". Figure |] shows the shapes of //s, fsD+^ fiNT+ and fiNT-- The 
SD^ component peaks at high muon and photon energy, making it the easiest of the SD 
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components to observe. This analysis, therefore, is mostly aimed at observing the SD"*" 
component. The form factors of the decay, Fy and Fa-, can, in principle, depend on g^, 
which is given by = Mj^- 2MkE^ in the K+ rest frame. In an 0{p^) ChPT calculation 
|PJ§], however, they are found to be independent and are given by Fy + Fa = 0.137, 
Fy — Fa = 0.052, which corresponds to BR{SD^) = 9.22 x 10^^. In the data analysis, we 
initially assume that they are constant, then test for dependence. 

The IB component of K^^^^ has been well measured in other experiments and found to 
agree with the QED prediction 0. The structure dependent components, on the other hand, 
have not yet been measured. For the SD^ component, the best limit is BR(SD''') < 3.0 x 10~^ 
0. There is also a limit on the combination BR(SD^+ INT^) < 1.3 x 10^'^ [0. In terms of 
the form factors, these limits translate into \Fy + Fa\ < 0.23, —0.3 < {Fy — Fa) < 2.5. 

Better measurements are available from the closely related process (Keuj)- 
In ^£1/7, the IB term is heavily suppressed by helicity, so that the SD terms are easier to 
extract. Since the INT terms are also highly suppressed, the signs of the form factors are 
practically impossible to measure. In O(p^) ChPT, Fy and Fa for K^wy are identical to 
those for K^^^. The K^^y experiments § give \Fv + Fa\ = 0.148 ± 0.010, \Fy - Fa\ < 0.49, 
in agreement with ChPT. 

In K^u^, the IB term is large, thus complicating the extraction of the SD terms, but 
also making the INT terms comparable in size to the SD terms. This makes it possible, 
in principle, to measure the sign as well as the magnitude of the form factors. In addition 
to its potential for checking the predictions of ChPT, K^^^ is also interesting as a probe of 
non-Standard Model CP- violation [^j. One can look for a T- violating component of muon 
polarization transverse to the plane of the decay. Such an effect is proportional to the INT 
components. 

The E787 experiment at the Brookhaven Alternating Gradient Synchrotron (AGS) [ |10|| , 
shown schematically in Figure |^, was used to look for the SD"*" component. E787, originally 
designed to search for —>-7['^uT7, uses a beam of mesons brought to a stop in a 
scintillating fiber target. From there, charged decay products can enter the drift chamber 
where their momentum is measured in a 1-T magnetic field. The charged tracks then enter 
the Range Stack (RS), which consists of 21 layers of scintillator and two layers of straw tube 
chambers (RSSC). Most tracks range out in the RS, thus allowing measurements of their 
total energy and range. A An photon detection system, composed of the Barrel Veto (BV) 
and two endcaps, surrounds the central region. In the present application the BV, covering 
70% of the solid angle and composed of lead and scintillator, is used to detect the photons 
of interest as well as to rule out the presence of more than one photon. 

The K^i,^ data was taken with the upgraded E787 detector, which was completed in 
1994. In this analysis, the redundant charged track energy and momentum measurements are 
combined (assuming a fi~^ mass) to give an improved measurement of the track kinematics. 
The rms resolution of this combined quantity is = 0.0164 ■ — 0.86MeV/c, for 205 < 
Pf^ < 236 MeV/c, where is the combined measurement expressed as a momentum. The 
resolutions for the azimuthal (0) and polar [6, with respect to the beam) angles of the muons 
are each 32 mrad. The resolutions on the photon kinematic quantities are cr^^ = 1.676^Je^ 
MeV {E^ in MeV), = 25 mrad, and erg = 45 mrad. 

A special trigger designed to search for the SD"*" component of K^^^ required a high 
energy charged track in the central region, a high energy photon in the BV, and no other 
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photons in the event. A two-day run using this trigger netted a total exposure of 9.2 x 10^ 
, yielding a total of 1.5 x 10^ K^^^ triggers. 

Analysis of the events passing the trigger proceeds in three steps: event reconstruction, 
background rejection, and K^^^ spectrum fitting. In the reconstruction step, the energy, 
time, and flight direction of the charged track and photon are calculated. Any additional 
photon energy not associated with the primary photon is also recorded. A kinematic fit to 
the K^j,^ hypothesis is applied to the charged track and the photon. Since there are four 
constraints (conservation of momentum and energy) and three unmeasured quantities (mo- 
mentum of the neutrino), the kinematics are over-constrained and Tion-K^^^ events should 
have a bad fit x^- Additionally, the kinematic fit yields measurements of and with 
better resolution than the raw quantities. These are the variables that are used in the final 
spectrum fits. 

The two main types of background that need to be rejected are K^2 accompanied by an 
accidental photon and i^"*" — ^vr^/i+z/^ {^^iz) or i^"*" — >■ tt+tt^ (ii',r2) where one of the photons 
from the 7r° decay satisfies the photon requirement and the other photon is undetected. The 
-f^^t2+accidental background can be suppressed in two independent ways: by requiring a 
tight time coincidence between the muon and the photon and by examining the kinematics 
of the decay. Since the accidental photon is randomly oriented relative to the muon, the 
cut on the of fhe kinematic fit to the K^j^^ hypothesis is especially effective against this 
background. Both K^^, and background can also be rejected in two independent ways: 
vetoing on any additional photon energy in the event and by the kinematics of the decay. 
The requirement that the charged track energy be above the K^-^ endpoint (-E'^+ > 137 MeV) 
is especially effective against this type of background. 

Since both types of background can be rejected by two independent methods, the total 
rejection of all cuts for each background may be calculated as the product of the rejections 
of the two methods. This allows a calculation of expected background based solely on the 
data, thus lowering the estimated systematic error. In the final signal region defined by 
E^+ > 137 MeV and E^ > 90 MeV, the expected background (with statistical error) from 
the i^'^2+ accidental source is 79.4 ± 4.8 events. The -K'^a and Kt^2 backgrounds are treated 
together and give a total expected background of 25.2 ± 3.8 events. 

Figure §(a) shows the final spectrum of events with the final signal region in the upper 
right corner delineated by the solid line. The number of events in this region is 2693, the vast 
majority of which are K^iy^. As a simple way of testing whether the K^iyj events are consis- 
tent with being only IB, we examine the distribution of the opening angle between the muon 
and the photon (cos6'^-y). Figure |^(b) shows this distribution for background-subtracted 
data. Superimposed on the data are Monte Carlo distributions for IB and SD^components 
of K^iy^. When only an IB component is allowed, the quality of the fit is very poor (x^= 
300, with 48 degrees of freedom). When an SD"*" contribution is allowed, a much better fit is 
obtained (x^= 58) [^, clearly indicating that a structure dependent component is present. 

The fit is incomplete, however, because it does not include the effects of the other 
fi^u^'-f components (SD~, INT"^, INT^). To include these effects, we generate Monte 
Carlo distributions with SD and INT components weighted by the form factors and normal- 
ized to the IB component. In Figure ^, we plot the between the E^+ vs. E^ histogram 
of this Monte Carlo sample and that observed in data (after background subtraction) as a 
function of the form factors. The minimum is 75 with 69 degrees of freedom. The best 
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fit values are 

\Fv + Fa\ = 0.165 ± 0.007, Fv-Fa^ 0.102 ± 0.073, 

where the errors are statistical. The minimum x^'s found in the regions where Fy + F4 < 
and where Fy + -F4 > differ by only 0.2. We thus have no information about the sign of 
Fy + Fa and can only measure its absolute value. The result corresponds to a branching 
ratio of BR(SD+) = (1.33 ± 0.12) x 10-^. 

The largest systematic errors associated with the form factor measurements come from 
possible distortions of the K^^^ spectrum induced by differences between the true detector 
and the Monte Carlo simulation. The two largest sources of distortion are non-linearity in 
the measurement of the photon energy and uncertainty in the thickness of the individual 
RS scintillator layers. For \Fy + Fa\-, these two sources lead to uncertainties of 0.0095 and 
0.0054, respectively. For \Fv — Fa\, they are 0.028 and 0.033. The systematic errors due 
to uncertainty in the level of background present in the final sample are estimated in the 
data-based background studies described above. They are found to be very small, totalling 
0.0007 for and 0.0097 for Fy — FA- Even a much enhanced background level would 

have only a small effect on the measurements. Adding the individual errors in quadrature, 
we find a total systematic error of 0.011 for \Fv -\- Fa\ and 0.044 for Fy — Fa 

As a check on possible systematic errors, the branching ratio for the IB component has 
also been extracted. This was accomplished by normalizing to a sample of K^2 decays that 
was taken simultaneously with the K^^^ data. For E"^ > 100 MeV and > 20MeV, we 
find BR(1B)= (3.6 ± 0.3) x 10~^, in good agreement with the theoretical value for this 
kinematic region, 3.3 x 10~^. Other checks included changing the binning of the vs. 
histogram and varying the cut. While both of these checks were limited by statistics, 
neither showed a systematic trend as the parameters were varied. Therefore, no systematic 
error is associated with these effects. 

As mentioned above, the form factors Fy and Fa have, to this point, been considered 
independent of . To assess the effect of including dependence, we assume the following 
form factor form: 

_ Fy{q^ = 0) _ FA{q^ = 0) 

^ 1 9 / 2 ' ^ 1 oil' 

L — q'^/my l~q^/m^ 

We take my = 0.870 GeV (the K* mass) and tua = 1.270 GeV (the Ki mass) and refit 
the measured ^^^1/7 spectrum in terms of the parameters Fy{Q) + Fa{0) and Fy{0) — Fa{0). 
The best fit parameters are: 

\Fy{Q) + Fa{Q)\ = 0.155 ± 0.008, Fy{0) - Fa{Q) = 0.062 ± 0.078. 

The corresponding SD"*" branching ratio is (1.37 ± 0.12) x 10~^. Although the value of 
|Fv'(0) + Fa{0)\ differs somewhat from that obtained assuming q'^ independence, the associ- 
ated branching ratio changes only slightly. Furthermore, the minimum of the fit is very 
insensitive to my and m^, so we are unable to measure them and cannot offer evidence of 
q^ dependence. 

In conclusion, we have observed a structure dependent component in the decay 
^ ^'^Vf/y. Under the assumption of (f independence, the associated form factors are 

\Fy + Fa| = 0.165 ± 0.007 ± 0.011, Fy - Fa^ 0.102 ± 0.073 ± 0.044. 
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Since the measurement of Fy — Fa is not significantly different from zero, we add statistical 
and systematic errors in quadrature and calculate the 90% confidence level: 



-0.04 < Fv- Fa< 0.24. 

The \Fv + Fa\ measurement is consistent with the previous result on — >■ but 
disagrees with the ChPT prediction by about two standard deviations. This is perhaps 

not surprising since at higher order in ChPT, kaon form factors are expected to differ from 
those of the pion [O. The 0{p^) calculation has been done for pions |]T3, but not yet for 



kaons. The limit on Fy — Fa is consistent with O(p^) ChPT and is significantly better than 
any previously obtained from kaon decay. A more detailed description of the analysis can 



be found in reference [0]. 
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FIGURES 




E, (MeV) (MeV) 

FIG. 1. Spectral shape functions for IB, SD+, INT+and INT" components of K^^^. The 
normahzations are arbitrary, and the scale on the fjB plot is logarithmic. The SD~ component 
is not shown because it peaks at low muon momentum and has negligible effect on the current 
analysis. 




FIG. 2. Side-view (a) and end-view (b) of upper half of the E787 detector. 
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FIG. 3. (a) Spectrum of events passing all but the final kinematic cut. At this point, nearly all 
events are i^^jy-y except those in the region < 137 MeV and > 150 MeV, where background 
from K^2, and i^,r2 is concentrated. The IB component of K^i,^ is concentrated at E^ < 100 MeV. 
The box marks the final cut of > 137 MeV and Ej > 90 MeV, within which the SD"^ component 
is enhanced, (b) Counts vs. cos(0ju-y) and various fits as described in text. 
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FIG. 4. contours for the fit to the E^+ vs. E~^ distribution, (a) Contours for all plausible 
values of the form factors. Each contour represents 50 units of x^- (b) Near a minimum. In 
this plot, each contoTir corresponds to one unit of x^- The one-standard-deviation uncertainties for 
Fy + Fa and Fy — Fa are also shown. 
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